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Rate of gravitational infiaton decay via gauge trace anomaly 
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We analyze decay processes of the infiaton field, <j>, during the coherent oscillation phase after 
inflation in f(4>)R gravity. It is inevitable that the infiaton decays gravitationally into gauge fields 
in the presence of f(<f>)R coupling. We show a concrete calculation of the rate that the infiaton 
field decays into a pair of gauge fields via the trace anomaly. Comparing this new decay channel 
via the anomaly with the channels from the tree-level analysis, we find that the branching ratio 
crucially depends on masses and the internal multiplicities (flavor quantum number) of decay product 
particles. While the infiaton decays exclusively into light fields, heavy fields still play a role in 
quantum loops. We argue that this process in principle allows us to constrain the effects of arbitrary 
heavy particles in the reheating. We also apply our analysis to Higgs inflation, and find that the 
gravitational decay rate would never exceed gauge interaction decay rates if quantum gravity is 
unimportant. 



I. INTRODUCTION 

Inflationary cosmology has passed a number of strin- 
gent observational tests, such as observations of cosmic 
microwave background temperature anisotropy [l|. De- 
spite the success of inflation models, the identity of the 
field that drives inflation is little understood. Any infla- 
tion model requires a graceful exit, the so-called reheat- 
ing after inflation (see reviews, e.g., [2 6] and references 
therein); otherwise the observed Universe cannot be pre- 
dicted by the models. The Universe must be thermalizcd 
and dominated by radiation fluids before the primordial 
nucleosynthesis. Thermalization is initiated by the de- 
cay of the infiaton field into lighter particles, and then 
the particles come to a state of local equilibrium with 
each other. 

It has been shown that inflation occurs naturally in 
models with nonminimal gravity [7H11J. and the spec 
trum of scalar curvature perturbations 
of tensor gravity wave perturbations [1 
fected by the presence of f(4>)R, thereby allowing us to 
constrain f(<j>) from the cosmological data. The f{4>)R 
gravity can also play an important role in the infiaton 
decay at the reheating |19l - l22l |. 

In theories with nonminimal couplings between the 
Ricci curvature and scalar fields, e.g., scalar-tensor grav- 
ity [H-H, Tt 2 gravity [HHH, supergravity [H-[3l|], and 
higher dimensional gravity theories, infiaton fields can 
decay via gravitational effects. The gravitational two- 
body decay rate is typically given by 
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where M Pl = (8irG)- 1/2 ~ 2.4 x 10 18 GcV, m a is the 
infiaton mass, and C is a mo del- dependent fudge factor. 
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The infiaton field condensate, a, decays into any pair of 
light fields if it is not conformally invariant. The un- 
known factor, C, significantly depends not only on the 
gravitational sector but also on the detailed properties 
of the matter sector, such as mass spectrum, spins, and 
the number of degrees of freedom at the energy scale of 
reheating. Thus, understanding physical grounds of C 
merits further understanding of the reheating. 

In this paper, we will take advantage of this gravita- 
tional nature in the reheating mechanism. We consider 
only the gravitational effect and assume no direct inter- 
action between the infiaton field and matter fields. Since 
gravity interacts universally, the infiaton field interacts 
with every field that exists at the reheating. Especially, 
we focus on an emergent interaction in the case where the 
infiaton couples to gravity nonminimally. Although the 
infiaton is assumed to be gauge-singlet, its decay prod- 
ucts (matter and radiation) may be charged under some 
gauge group, e.g., SU(N). We aim to constrain phys- 
ical properties of new particle species at the reheating, 
such as spin, mass spectrum, and the number of degrees 
of freedom. As was discussed in (20[, the large number 
of degrees of freedom significantly enhances the gravita- 
tional decay. 

Our approach can be contrasted with the recent devel- 
opment of the Higgs inflation model [32M361 ] . where no 
new degrees of freedom are introduced. The model with 
nonminimal kinetic gravitational coupling, G^d^d^c/), 
is another possibility [37, 38] in this direction. Reheating 
in Higgs inflation has been considered by [H, E3l • 

The organization of this paper is as follows. In Sec. |TT] 
an inflation model with nonminimal gravitational cou- 
pling and t/(l)-charged matter Lagrangian is given. The 
gravitationally induced interactions are derived. The 
rates of gravitational infiaton decay are given in Sec. IIIII 
Our main result, the decay rate via conformal (trace) 
anomaly [Eq. (|25p] , is presented in subsection IIII B[ Im- 
plications of the gravitational infiaton decay for reheat- 
ing are investigated in Sec. IIVI An implication for Higgs 
inflation is given there. We show a detailed deriva- 
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tion of fermionic one-loop decay rate via trace anomaly 
in the Appendix. We work with the metric signature 
(+, —,—,—) throughout the paper. This sig n convention 
is different from that in previous works |l9l |20| . 



II. f(<p)R GRAVITY AND INDUCED 
COUPLINGS 



For concreteness, we assume that the inflaton, <fi, is a 
neutral scalar field that is nonminimally coupled to grav- 
ity but gauge-singlet, i.e., no direct interaction with the 
SU(3)c x SU(2)l x U(l)y standard model sector or any 
other matter fields. To make physics clear, we employ a 
QED-like model, [/(l)-charged scalars and fermions, as 
a matter sector at the energy scale of reheating: 
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where N x and are the internal flavor quantum num- 



ber for scalar and fcrmion, respectively. We impose the 
boundary condition f(v) = A/p ; , where v is the vacuum 
expectation value (vev) of <fr at the potential minimum, to 
guarantee the ordinary Einstein gravity at low energy. 1 
We have explicitly shown the square root of the determi- 
nant of the metric, \J—g, to emphasize that two sectors 
talk to each other only minimally via gravity. The co- 
variant derivatives for each field arc defined as 



D» x = g^{d v -igA v ) X , 
T^> = e" a 7 a (a„-r M -iflA M )V, 
F» v = g^g va F ptT = g^g va (V p A a - V a A p ) 
= g w g va {d p A a -d a A p ), 



(3) 



where g is the gauge coupling constant and assumed to 
be weakly coupled, g < 0(1). Here a, f3, 7, • • • denote 
Lorentz indices while A, /j,, u, ■ ■ ■ denote general coordi- 
nate indices, a is a tetrad (vierbein) field. Y ^ is a spin 
connection and S"^are generators of the Lorentz group 
-iE Q/3 e A Q V A1 e A/3 and E Qf3 = -YP a = 
e.g., Sec. 12.5 in g2; Sec. 3.8 in g|). 
This simple model is rich enough to demonstrate impor- 
tant gravitational decay processes in reheating as we will 
see in the following section. 



given by 1^ e 
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1 This dynamic boundary condition is motivated by the spirit of 
Zee's induced gravity theory [4lj | and also applied to other models 
of nonminimally coupled inflation 0, l2lf . 



After inflation the inflaton field oscillates about the 
minimum of the potential. Thus we expand 4> as 



v + a, 



where a represents coherent condensate of inflaton 
quanta measured from its vev, (</)) = v. In the Jordan 
frame, where the theory was originally defined, the in- 
flaton, a, decays into the matter sector through loops 
and mixing involving graviton if there is no direct in- 
teraction. It is simpler to compute physical quantities 
in the Einstein frame, where inflaton and graviton are 
diagonalized. 

Performing conformal transformation to the metric 
tensor 
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f(d>) 
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and rescaling fields simultaneously, 2 3 

x -> x = n~ 1 x, 



(T — > a = <j\ \ 1 



3 ( f'(v) 
2 I Mpi 



ij) -> i> = rr 3/2 v>, 



A^ -> Afj, = A^ , A" -> A" = Q- 2 A», 



(6) 



interactions are spontaneously induced by gravity. 
Transformed matter Lagrangians are given by 



4 = W^M'Tv**' 



xtu, (7) 
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where the covariant derivative for scalars is defined as 

V^x^d^x + xd^lnty-igA^x- (10) 



2 Since canonical kinetic terms are needed to quantize fields in the 
new coordinate (Einstein frame), the field redefinitions JBJ are 
necessary. The rescaling of a follows immediately from l44f 



4>^<t> = M Pl 



1 3 f/W 2 
/(*) 



3 In the linearized theory, infinitesimal conformal transformation 
simply corresponds to diagonalization between graviton and in- 
flaton, cr, canonically normalized. Since two frames are very 
close at the potential minimum of the inflaton, we do not have 
to rescale matter fields [TJJ. If one would like to know nonlin- 
ear interaction, finite conformal transformation is necessary and 
matter fields should also be rescaled accordingly. 
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Prom Eqs. 0, © and © it can be seen that masslcss 
fermions and gauge bosons are conformally invariant at 
the classical level while massless scalars are not. 4 

Interaction between the inflaton field and matter fields 
appears through the conformal factor Q(<t), 



X+ 



f'(v)a , f"(v)a 2 , f"'(v)* 3 
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m Pl 
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where we have defined 



6M p; 



,(H) 



Fn(«) = 



d n f/d(j> n (v) 



[f'(v)/M, 



pi 



T./2 ' 



(12) 



where n is the integer number and f^ n \v) = d n f/d4> n (v) 
is the n th derivative of f((f>) with respect to <f> at vev. 
Since a/Mpi <C 1 during the coherent oscillation phase 
after inflation, it is natural to assume the series (jTTj) is 
convergent. 

The interaction Lagrangians to the lowest order in a 
are given by expanding Eqs. © and ©: 



pi 

+9> iv {d tl &)tD v x + 2m 2 x ax*x) 
^^Mf l i^d^irx) + 2rn 2 x ax*x) (13) 

^HTFT^ H""^*)* D vX + 2mlx*x) (14) 
iv1 pi 



r - 



9 KT72 <rW- 



2M 



(15) 



pi 



In deriving the last line of C axx , we have used integration 
by parts and the classical equations of motion for x an d 



4 The spinor and gauge field connections are conformally invari- 
ant: 

fV = -iE a ' 8 e A a V M (e A/3 ) 



i 



Ffiv = Ffiu, F^ = n-^F^. 



G oooooooooooo 



tj oooooooooooo 




FIG. 1: Gravitational inflaton decay at tree level. 



X* from Eq. ©: 



MV-gg^d^x) 



-g ( r 2 raj + g 2 A 2 - 2igA^d v ) X, 



duW-gg^d^x*) 



g [ n- 2 m 2 + g 2 A 2 + 2igA^d v ) x* , (16) 



where d^A 11 = has been used. The classical field equa- 
tions correspond to the on-mass shell condition in the 
Feynman diagrams. Therefore, the expression (|14j) can 
be used only for analysis at tree level (see Figs. [I] [2]). F° r 
an analysis off the mass shell one has to use the expres- 
sion (fl3|) instead (see Fig. [3]). Note that neither 4- leg 
nor 5-leg interaction appears due to exact cancellation 
in the a(D M x)* D^X term on the mass shell (see Fig. @|. 
Note also that higher order interaction Lagrangians can 
be derived systematically by expanding Eqs. © and © 
in power scries of a. In the following sections we com- 
pute exclusively in the Einstein frame. We shall suppress 
carets on variables when it does not cause any confusion. 



III. INFLATON DECAY BY BREAKING 
CONFORMAL INVARIANCE 



A. Induced decay at tree level 



With the interaction Lagrangians, C axx and C a ^ !1 p 1 the 
inflaton decays into a light pair of scalers or fermions at 
tree level (Fig. [T]). The rates are given by the standard 
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FIG. 2: Three-body decays at tree level. 



quantum field theory analysis [41 

rO -> x+x-) 



1 



4m, 



1/2 



(17) 
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, (18) 



where gravitationally induced coupling constants are de- 
fined as 



.'A 

<7y 



Fi(t>)(m£ + 2m^ 



4M* ; 



_ Fi(v)m v 



2M% 1 



(19) 
(20) 



Note that g x is evaluated on the mass shell, and the first 
term is due to the derivative coupling of scalars. 5 

The tree-level two-body decays are kinematically sup- 
pressed as seen in Eqs. (fTTl) and ([T8|) . Especially, the 
inflaton cannot decay into scalars or fermions if decay 
products are heavier than half mass of the inflaton, 
m a < 2m x , 2m w . It is usual to have the global fla- 
vor symmetry broken, and masses of particles arc differ- 
ent as in particle physics. Thus, we split flavor quan- 



tum number into two pieces: N x 



N. 



N, 



N Xh and 



N^ h , where N Xl and N wt denote the num- 



ber of species lighter than half of the inflaton mass while 
N Xh and N^ h denote the number of heavier species. In 
this case one should replace N x , in Eqs. (fTT|> and (JTT 



with N Xe , N^ e 



At tree level the inflaton also decays into three-body 
(Fig. [2]) and four-body final states. Gauge bosons are pro- 
duced in the radiative decay (bremsstrahlung) processes. 
From Eq. ([14"]) one may naively expect scalar QED-like 
interactions shown in Fig. |4j However, they do not ap- 
pear if one evaluates the derivative couplings properly. 

How fast do these processes proceed? Three-body de- 
cays are phase-space suppressed compared to the two- 
body decays. For bosonic decays 



IV x+x-) 



N x g x 

47rm (T 
a 



8tt 



64?rM^ 



r(o- x+x-), 



(21) 



where m x <C m 
fermionic decays 



r(cr -> iplp) 



is assumed and a = g 2 /(47r). For 



8?r 
— 1 (er 



32ttM£, 



(22) 



where mw, <C m CT is assumed. If g ~ O(0. 1), r(cr — > 
3)/r(<7 -> 2) ~ 5 2 /(32tt 2 ) ~ 0(1O'- 4 - 10- 3 ) for both 
X and ?/>. 6 Although the branching ratio of three-body 
decay is small, pairs of decay products quickly annihilate 
into gauge bosons right after the gravitational inflaton 
decay. The annihilation process proceeds much faster 
than the gravitational decay if g ~ O(0. 1). 

In four-body decays, the inflaton radiates two gauge 
fields or another pair of charged fields in addition to a 
pair of charged fields. We have assumed that the gauge 
interaction is the only renormalizable interaction of the 
matter fields, x an d V'j n0 Yukawa interaction and no 
self-interaction arc introduced. If there were the Yukawa 
interaction in the original frame, yx^ip + yx*'4 ) '4 ! i then 
the four-leg interaction, y&x^tp + y&X*' , l J '4>> would show 
up in the Einstein frame. If there were the scalar self- 
interaction, A|x| 4 , the five- leg interaction, Act|x| 4 , would 
show up. As is the case of three-body decays, four-body 
decays are phase-space suppressed compared to two-body 
and three-body decays. 



B. Induced decay via quantum trace anomaly 

In this subsection, we present the main result of this 
paper: the rate of gravitational inflaton decay into gauge 
fields via gauge trace anomaly [Eq. ([2"5]l]. 

At the classical level the inflaton can decay into only 
light degrees of freedom. At the quantum level, however, 



5 For the scalaron decay in Starobinsky's 7?. 2 -innation [231428H , one 
can identify F±(v)/Mpi = 2/v / 6 and m a = [i for the scalaron 
mass. 



° Soft photon processes in Fig. [2] are known to be infrared log 
divergent. Therefore, we consider only hard photon processes in 
Fig. [2] by introducing a finite infrared energy cutoff for outgoing 
photons. 
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the inflaton can also decay into heavy degrees of freedom 
if they are intermediate states. In Fig. [3] we show the 
two gauge boson decays of the inflaton via fermionic and 
bosonic one-loops. Even though there is no induced cou- 
pling between the inflaton and massless gauge bosons, 
the inflaton still decays into gauge bosons through the 
anomaly process. 

In the Einstein frame the fermionic matter field, ip, 
couples to the inflaton by g^ = Fi(v)m^ / (2Mp Y ) , and it 
also couples to gauge fields minimally. The fermions me- 
diate between a and fields as a triangle loop (Fig. [3]). 
We find that the decay rate via fermionic one-loop is 
given by 



r f {v->2A tl ) 



256n 3 M^ 



(23) 



oooooooooo 
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r 



A v (k') 




A„{k') 



where a = g 2 /(47r). It is remarkable that the intermedi- 
ate particle mass dependence appears only in the func- 
tion 1 f (m 2 1 'm 2 ^) . A detailed derivation is given in the 
Appendix. If the ip field is heavier than the inflaton, the 
tree-level decay [Eq. (fTS]) ] is kinematically forbidden. But 
the decay channel to gauge bosons is still open through 
heavy fermions. 

The decay rate via bosonic one- loop (Fig. [3]) is similarly 
given by 



T s {a^2A^) 

10247t 3 M4j 



.(24) 



The functions, If and / s , are given in the Ap- 
pendix. In the heavy intermediate particle limit we 
have Ifiml/m 2 ^ -> 0) = 1/3 and l s {m 2 a lm 2 x -> 
0) = 1/6, and in the light intermediate particle limit 
//(m^/m^, — s- oo ) = l s {m 2 a lm 2 x — s- oo) = while 
(2 + m 2 r /m 2 ,)I s (m 2 r /m 2 ,) —> 2. Thus, only heavy interme- 
diate particles contribute to the fermion loop [Eq. (|23|) ] 
while both light and heavy ones do to the scalar loop 
[Eq. (23)]. 

The total decay rate for the process a — > 2A ll is in 



general given by 



r(<7 -> 2A^) 



10247r 3 M^ 




+ 



E 

s=l 



(25) 



We have assumed that there is no massive gauge field at 
the reheating. If the energy scale of reheating is lower 
than that of the electroweak scale, three of four massless 
and Wt bosons become massive and bosons 
since the SU(2) L x U(l)y gauge group is spontaneously 
broken down to U(1)em- To include a massive gauge 
boson in our model is straightforward but beyond the 
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FIG. 3: Gravitational inflaton decay at one-loop. 



scope of this paper. The vector boson loop contribution 
might be as significant as scalars and fermions, but its 
evaluation is more involved than the present model since 
one has to include the Faddeev-Popov ghost field in the 
loop diagrams [46|, |47j . We will briefly discuss an effect 
from massive gauge bosons in the following. 

What is the physical origin of these decay processes? 
If there is a dimensionful parameter, such as mass, in the 
Lagrangian, it breaks scale invariance explicitly. Confor- 
mal invariance is also broken because the scale symmetry 
is a subgroup of the conformal symmetry. The emergent 
interactions [Eqs. (TT3"]) and (fI5]l] can be understood as 
a result of breaking this conformal symmetry. At the 
classical level there is no (energy) scale dependence in 
the dimensionlcss parameter, such as the gauge coupling 
constant. At the quantum level, however, the conformal 
invariance ought to be broken as the gauge coupling con- 
stant, g, runs with its energy scale. 

The breaking of conformal invariance leads to nonzero 
value in the trace of the classical energy-momentum ten- 
sor El: 



rpp, 

J m(" 



n ss n 



(26) 
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In deriving this expression, we use definitions of the 
energy-momentum tensor and functional differentiation 
with respect to Q: 



— S m + 

- s m + 

where we have used 

rpm 
{11/ 

bp* 



Fi(u)5- 



(27) 



2 5S m 
f=g Sgv ' 
-2g» v 6lnfl, 



(28) 



in the second line of Eq. (|27j) . and have used Eq. (fTTj) 
in the last line of Eq. (f27|) . The inflaton field, tr, thus, 
couples to the matter fields via the trace of the energy- 
momentum tensor of matter fields. Suppressing carets 
on variables again, the leading effective interaction La- 
grangian can be written as (cf. (4(| [4?J ) 

C . t = V^^T- 
2M 2 



pi 



= Y j n-(D»Xs)*D» Xs + 2U(x* sX s)} 



/ = 1 



v=l 



2g 



F^F^, (29) 



where indices /, s, and v run through massive fermion, 
scalar, and vector field species, respectively. We have in- 
cluded the contribution from massive non-Abelian gauge 
fields. Here fth(g) is the heavy particle contribution of the 
Callan-Symanzik (or Gell-Mann-Low) ft function and is 
given by (cf. (47j ) 



h(g) 
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(4tt) : 



E 

heavy 









I) 







1 



-N f 



+0(g% (30) 



where N s , Nf, and N v are the internal (flavor) quantum 
number of scalar, fermion, and SU (N)-chaxged massive 
vector species, respectively. The term, —N v /3, inside the 
bracket stands for longitudinal components of massive 
vector fields. The summation is taken over all species 
heavier than the inflaton, and their masses should be less 
than a cutoff scale of the model. This amounts to only 
the one-loop correction to the effective action. The higher 
order loops can be ignored if the perturbation expansion 
converges: g < 1. Since the running of g is determined by 
effectively massless particles, fth(g) does not contribute 
to the running of g until the energy scale reaches the 
ultraviolet region, \x > rriheavy 3> 



Equation (|29|) takes into account only U(l) charge 
renormalization as the leading effect to produce pairs 
of massless U(l) gauge bosons. For a consistent one- 
loop analysis, one needs to include charge renormaliza- 
tion of SU(N) gauge fields, ft(g v )V^V*^ /2g v , where 
V£, = W - dvV£ +9vfabcV»V„ c and the indices a, 6, c 
take the value +, — , 3, • ■ • ,N 2 — 1. Renormalization of 
massive vector fields is beyond the scope of this paper. 

In the heavy intermediate particle limit, nii 3> m a , the 
decay width can be computed from the effective interac- 
tion Lagrangian [Eq. (|2T)|) ]. The inflaton decays into a 
pair of massless gauge fields with 



2A„ 



1024tt 3 MJ, ; 



i=s,f,v 



(31) 



where bi are the first coefficients of the ft function [inside 
the bracket of Eq. (f3T))) ] . The summation is taken over all 
species heavier than the inflaton. With no massive vector 
species, N v = 0, the expression agrees with Eq. (|2l)|) in 
the limit of rrii 3> m a . 



IV. IMPLICATION FOR THE PHYSICS OF 
REHEATING 



What could these decay rates imply for the physics of 
reheating after inflation? In general, the gaugc-singlct- 
inflaton field has difficulty in reheating. However, the 
presence of the f(<f>)R term enables the singlet-inflaton 
to interact with matter as was shown in the previous 
sections. Therefore, the Universe reheats naturally [lj|. 
In order to connect the gravitational inflaton decay to 
the observable Universe, one actually needs to know the 
matter contents of the Universe after inflation. Some of 
them are associated with a visible sector involving the 
standard model particles. Others are associated with a 
hidden sector involving dark matter. While we have not 
addressed these distinctions, new hints are provided by 
our model; the gravitational decay rate crucially depends 
on spins, mass spectra, and the number of degrees of 
freedom. 

As the rates depend on the number of species in a spe- 
cific way, they in principle allow us to constrain the num- 
ber of heavy particles at the reheating [Eqs. (| 1 T[) and (fT5)) 
for light particles, Eqs. (|2"3"]l and ([M| for heavy particles]. 
Let us split the degrees of freedom into light (me -C m a ) 
and heavy (rrih m a ) species as A^ x 



N xe + N Xh and 



-ZV,/, = N^p t + N$ b by assuming flavor symmetry broken 
into two pieces. For scalars, in the absence of fcrmions, 
the loop process [Eq. (|24[) ] becomes faster than the tree- 
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level process [Eq. (|XT[)] if 

r s (<7^2^ M ) _ a 2 (6A x ,+A x J 2 



r(<7 -> x+x-) 

1 /TV 
> 120- ^ 
Xh ~ a V 10 



144tt 2 A y 



> l; 



1/2 



if A Xh »A Xf , 



^^40— if A Xh «A x ,. 



(32) 

(33) 
(34) 



Two remarks on the condition must be made. The first is 
that the channel to light scalars, including the standard 
model Higgs bosons, dominates over the anomaly channel 
if the condition [Eq. ([32]) ] is not met. The second is that 
Eq. (|34[) holds only if the gauge interaction remains weak; 
it might hit the strong coupling scale with sufficiently 
large N XI below or at the reheating energy scale, T r h, 
since all effectively massless charged scalars contribute to 
the running of gauge coupling with the positive sign. Let 
us assume the hierarchy of scales as fi ~ T r h <C m a . (If 
T r h 3> rricr, the thermal effects are important (l9j.) Then, 
if me <C T r h, N Xe species contribute to ft function and the 
gauge coupling would become strong at T r h spoiling the 
perturbative analysis made in this paper. If mi 3> T r h, 
N xe species do not contribute to (3 function and the gauge 
coupling remains weak at T r h- The presence of heavy 
states does not matter in driving to the strong coupling 
scale but in enhancing the anomalous decay process. 

Similarly for fermions, in the absence of scalars Xsi the 
loop process [Eq. ([2"3"jl ] becomes faster than the tree-level 
process [Eq. ljlg]t ] if 



12-K 2 m\N^ 



> 270- 



100 



> 1; 



1/2 



(35) 
(36) 



which can be satisfied more easily than bosonic processes 
since m 2 ^> to 2 ^ . Therefore, the large number of heavy 
fermions may affect the reheating process significantly. 

In the presence of both scalars and fermions, by com- 
paring Eq. (f25j) to Eqs. (flTf and (fl8]l. one finds the con- 
dition on which the anomaly induced decay is important: 



r(cr -> 2A^) a 2 (6N X 



7V V 



r. 



1447T 2 



N XI + 2^fN € 



(37) 



If there are many charged heavy species at the reheat- 
ing, gauge fields can be generated efficiently via trace 
anomaly. Moreover if the reheating process is dominated 
by this channel, one can constrain the number of heavy 
species by the reheating temperature as in the argument 
of [lj|; too high reheating temperature would produce 
unwanted relics, e.g., topological defects in the grand uni- 
fied theory. 

Finally, let us mention the gravitational decay in the 
Higgs inflation model [32[. The Higgs inflaton directly 



couples to the matter with the gauge and Yukawa in- 
teractions. Since the vev of the Higgs field is small 
compared to the Planck scale in any sensible gauge 
theories, the gravitationally induced couplings are set 
to be small: 7 g x ~ £i>m 2 /(4Mp ; ) -C gm x /mv and 
~ £,vm^/(2Mp l ) <C y ~ gm^/my < 0(1), where 
m a and my are masses of the Higgs and massive gauge 
boson in the Einstein frame, respectively. Also, g and 
y denote gauge and Yukawa coupling strengths, and are 
not much smaller than 0(1). For simplicity, we compare 
perturbative decay rates for fermions: r*™° /r*™v ~ 
y 2 l ML/{^ 2 v 2 m tT m^ t ) ~ g 2 m^ t M%J(^v 2 m a m\) > 
0(1). 8 If the process is kinematically forbidden, the one- 
loop effect takes place instead: r g ^- loo P/r g ™- loop ~ 
VfcAfpj/(£Vm^) ~ S »4/({V<)» 0(1). In both 
cases the ratios exhibit the hierarchy between gauge force 
and gravity. The main energy transfer, therefore, un- 
likely goes through gravitational decays unless quantum 
gravity becomes important. Note that the cutoff scale 
is given by A ~ Mpi/{\N^), where N is the number of 
species. 

So far the analysis has been limited to renormalizable 
operators and a few of the dimcnsion-5 operators, such 
as tj,va(dx) 2 and £,vaF 2 . The smallness of v made those 
decay channels inefficient in the standard model. In fact, 
the Higgs inflation model contains up to dimension-6 op- 
erators induced by nonminimal gravitational coupling, 
and we have not exhausted all of them. The operators, 
such as £cr 2 x 2 , ^ct 2 ^^, £,& 2 (dx) 2 , £,& 2 F 2 , do not include 
small Higgs vev in the induced couplings, and thus might 
dominate over the decays from operators considered in 
this work. The gravitational pair annihilation rate of in- 
flatons to scalars by the dimension-6 operators was cal- 
culated in 1201. 



V. CONCLUSION 

We have studied the gravitational inflaton decays in 
f(4>)R gravity with the renormalizable QED-like matter 
sector. The model does not require any direct couplings 
between the inflaton and matter fields. We have shown 
that the inflaton must decay into massless gauge fields via 
the trace anomaly process with the rate of Eq. (|25p. The 
decay channel is interesting because its amplitude is sen- 
sitive to the contribution of all the charged species in the 
theory including very massive ones. Especially if the the- 
ory includes the large number of heavy charged species, 
they do not appear as main compositions of the primor- 
dial radiation plasma but play the role of the shadow Cab- 



frcc 



7 On the contrary, the vev of the gauge-singlct-infiaton is 
parameter, which can be ~ Mpi . 

8 In the standard model, in fact, the decay proceeds nonperturba- 
tively due to violation of the adiabaticity condition, and is more 



efficient than perturbative one 39, 4i 
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FIG. 4: Naive diagrams at tree-level are forbidden. 



inet enhancing the gravitational decay into gauge fields 
enormously. 

Since the gravitational couplings are induced by break- 
ing of conformal invariance, it affects all of the nonconfor- 
mal fields that exist at the reheating. We have considered 
implications of the anomaly induced decay for the physics 
of reheating in a simple case, where the matter particles 
are either light (m, <§; m a ) or heavy (nii m a ) com- 
pared to the inflaton. The fermion loop process depends 
only on the number of heavy species while the scalar loop 
process depends on that of both light and heavy species. 
This new decay channel may rescue some inflation mod- 
els that fail to reheat sufficiently. If an inflation model 
requires many degrees of freedom beyond the standard 



model, the gravitational decay via trace anomaly would 
contribute substantially. Our argument would help to 
constrain the nonminimal particle physics of reheating. 
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Appendix A: Calculation of T-matrix amplitude for fermionic one-loop 



We use the dimensional rcgularization scheme, and then the amplitude of the fermionic triangle diagram is finite. 
The transition matrix (T-matrix) amplitude for a — > 2A^ via the fermionic one-loop is 

(A(k)A(k')\T\a(q)) = {2n) i S 4 (q- k- k')M, 

, • ., . 2 Ar *, I v.,,,v/" d 4 s Tr[(^+m) 7 ^+ ^ + m) 7 I/ (./+ d+m)} 

M = Mx+M 2 = ig 2 9 f N f e;{k)e:(k') / 7^7^ / f . v /'^ 2 , • u/ , ^2 2 , m + fc fc 

p J (27r) 4 (s z — m z + ie)((s + k) z — m z + ie){{s + q) z — m z + ie) 

= lg z g f N f [e;(k)et(k')I^(k 7 k') + e;(k')et(k)I^ (k' , k)], 



where we have used the Feynman slash notation for 4- momentum, <f = 7 M g M . The second term in Ai represents the 
triangle diagram with the opposite charge current (or internal momentum) direction, <?/ is a gravitationally induced 
coupling constant between a and an intermediate fermion particle, and m = mj is the mass of the intermediate 
particle. 

To carry out loop momentum integration we need to combine propagators. Using the Feynman parameter trick, 



9 



we rewrite the denominator as 
1 



D (s 2 — m 2 + ie)((s + k) 2 — m 2 + ie)((s + q) 2 — m 2 + ie) 

dx I dyx{(l — x)(s — m + ie) + xy((s + k) —m + ie) + x(l — y)((s + q) — m + ie)} 



where 



dx dyx{£ 2 - A 2 + ie}-'\ (A2) 



i" = + k> i x + k' f *x(l-y), A 2 =m 2 -fx 2 x(l-x)(l-y). (A3) 



Here we have used q^ = k^ + k ,fJ " and on-mass shell conditions: k 2 = k' 2 = 0, q 2 = fi 2 (p 2 = m 2 , in the rest frame of 
the inflaton). The numerator in the loop integral is 

N^ v = Tr( 4 + m)^ (4 + # + m)^ (4+4+ m) 

= Am^s" + 2{s^k v + k^s u + s^q v ) + k^q v + q^k v 

+(-2s-k-k-q- s 2 + m 2 )g^], (A4) 

where we have used the 7-matrix algebra 

Tr(ocM#7) = 0, 

Tr^^V) = K9^9 af> ~ S^V" + 9^9™)- (A5) 



Shifting the loop momentum by Eq. (|A3I) . one can rewrite the numerator as 

= Am[A£^t + A{k^x + k' tl x{l~y)){k' J x + k lU x{l~y)) 

-2x(3k^k v + Wk'") - 2x{\ - y)(2fc"V + k^k' v + k llL k' v ) 
+2k"k v + k^k' u + k^k v 

+ {2x(l - y)pi 2 - fj, 2 /2 -£ 2 -x 2 (l- y)v 2 + m 2 )g^} 
= 4m[2(l - x)(l - 2x)k^k v - 2x{\ - y)(l - 2x{l - y))k""k" J 
+ (1 - 2a:) (1 - 2x(l - y))k^k' v 
+ {l-Ax{l-x){l-y))k^k u 

-(1 - 4a:(l - x){\ -y) + 2x 2 {\ - y) VV72 + mVl, (A6) 
where the expression is valid inside the loop integral, and we have used 4-dimcnsional integrals in Minkowski space 

MH" f d 4 £ t 2 g» v 



/ 



(2tt) 4 {i 2 - A 2 + ie} 3 J (2tt) 4 {I 2 - A 2 + ie} 3 
d 4 £ £^ 



(2tt) 4 {£ 2 -A 2 + ie} 3 



0. (A7) 



Since the integrals converge, one can set d = 4 from the start. Thus, we have evaluated the 7 matrices and loop 
integrations in 4-dimension. 

The invariant matrix amplitude must be gauge invariant in a given order of perturbation theory. The Ward identity 
implies 

k^M^ = k' u M^ = 0, 

fc,Mr = -VWf, k'„MT = -KM 1 ?. (A8) 
The invariant matrix amplitude, therefore, becomes 

Mi = ig 2 <j f N, e ;,(k,\K(k',\')ir, 



10 



where we have explicitly shown polarization of gauge fields, A, and have used 

d 4 £ 1 -i 



(2tt) 4 {£ 2 - A 2 + ie} 3 2(4tt) 2 A 2 



(A10) 



A function, Ii([i 2 /raj), represents mass or energy dependence of an intermediate particle. For fermionic one-loop, 



dx I dyx 



dx 



dy 



1-42(1-30(1-2/) 

l-£c(l-s)(l-y) 
1 — 4xy 
1 - £xy 



arcsm 



2 VI 



1-i 1- 



In 



l-Vl-4/£ 





2" 


j — 27T 





if £ < 4 

if e>4 



This result agrees with [46H50J ■ The closed form formulas are first presented in |47| . 
Similarly for scalar one- loop (cf. [4ol |47j). 



Is(0 



dx 



l-x 



dy 



Axy 
1 - £xy 

_ |arcsin 2 (^) 
In 



if £<4 



The decay rate of a — > 2A I1 can be calculated as 

1 1 1 



if C>4 



(2tt) 2 2/i 2 



64tt 3 



at, 



^- — * 777, j 



(All) 



(A12) 



(A13) 



where a = g 2 /(Air). Note that in the heavy intermediate particle limit, we have 7/(£ — s- 0) = 1/3 and I s (^ 0) = 1/6, 
and in the light intermediate particle limit, — > 00) = J s (£ — > 00) = 0. One can evaluate massive vector loop 
diagrams similarly; the rate should be of the same order of magnitude as the fermion and scalar loop diagrams. 



[1] E. Komatsu et al. (2010), 1001.4538. 

[2] A. D. Linde, Particle Physics and Inflationary Cosmology 
(Harwood, Chur, Switzerland, 1990), hep-th/0503203. 

[3] D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999), 
hep-ph/9807278. 

[4] V. Mukhanov, Physical Foundations of Cosmology (Cam- 
bridge University Press, NY, 2005). 

[5] S. Weinberg, Cosmology (Oxford University Press, NY, 
2008). 

[6] J. Braden, L. Kofman, and N. Barnaby, JCAP 1007, 016 

(2010), 1005.2196. 
[7] F. S. Accetta, D. J. Zoller, and M. S. Turner, Phys. Rev. 

D31, 3046 (1985). 
[8] D. La and P. J. Steinhardt, Phys. Rev. Lett. 62, 376 

(1989). 

[9] T. Futamase and K.-i. Maeda, Phys. Rev. D39, 399 
(1989). 



[10] D. S. Salopek, J. R. Bond, and J. M. Bardeen, Phys. Rev. 

D40, 1753 (1989). 
[11] R. Fakir and W. G. Unruh, Phys. Rev. D41, 1783 (1990). 
[12] N. Makino and M. Sasaki, Prog. Theor. Phys. 86, 103 

(1991). 

[13] D. I. Kaiser, Phys. Rev. D52, 4295 (1995), astro- 
ph/9408044. 

[14] V. Faraoni, Phys. Lett. A269, 209 (2000), gr- 
qc/0004007. 

[15] E. Komatsu and T. Futamase, Phys. Rev. D58, 023004 

(1998) , astro-ph/9711340. 

[16] E. Komatsu and T. Futamase, Phys. Rev. D59, 064029 

(1999) , astro-ph/9901127. 

[17] J. C. Hwang and H. Noh, Phys. Rev. Lett. 81, 5274 
(1998). 

[18] S. Tsujikawa and B. Gumjudpai, Phys. Rev. D69, 123523 
(2004), astro-ph/0402185. 



11 



[19] Y. Watanabe and E. Komatsu, Phys. Rev. D75, 061301 

(2007) , gr-qc/0612120. 

[20] Y. Watanabe and E. Komatsu, Phys. Rev. D77, 043514 

(2008) , 0711.3442. 

[21] N. Kaloper, L. Sorbo, and J. Yokoyama, Phys. Rev. D78, 

043527 (2008), 0803.3809. 
[22] C. Pallis, Phys. Lett. B692, 287 (2010), 1002.4765. 
[23] A. A. Starobinsky, Phys. Lett. B91, 99 (1980). 
[24] A. Vilenkin, Phys. Rev. D32, 2511 (1985). 
[25] M. B. Mijic, M. S. Morris, and W.-M. Suen, Phys. Rev. 

D34, 2934 (1986). 
[26] S. Kalara, N. Kaloper, and K. A. Olive, Nucl. Phys. 

B341, 252 (1990). 
[27] T. Faulkner, M. Tegmark, E. F. Bunn, and Y. Mao, Phys. 

Rev. D76, 063505 (2007), astro-ph/0612569. 
[28] D. S. Gorbunov and A. G. Panin (2010), 1009.2448. 
[29] M. Endo, M. Kawasaki, F. Takahashi, and T. T. 

Yanagida, Phys. Lett. B642, 518 (2006), hep- 

ph/0607170. 

[30] M. Endo, F. Takahashi, and T. T. Yanagida, Phys. Rev. 

D76, 083509 (2007), 0706.0986. 
[31] M. Endo, F. Takahashi, and T. T. Yanagida, Phys. Lett. 

B658, 236 (2008), hep-ph/0701042. 
[32] F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B659, 

703 (2008), 0710.3755. 
[33] A. De Simone, M. P. Hertzberg, and F. Wilczek, Phys. 

Lett. B678, 1 (2009), 0812.4946. 
[34] F. Bezrukov and M. Shaposhnikov, JHEP 07, 089 (2009), 

0904.1537. 

[35] A. O. Barvinsky, A. Y. Kamenshchik, C. Kiefer, A. A. 



Starobinsky, and C. F. Steinwachs (2009), 0910.1041. 
[36] F. Bezrukov, A. Magnin, M. Shaposhnikov, and 

S. Sibiryakov (2010), 1008.5157. 
[37] C. Germani and A. Kehagias, Phys. Rev. Lett. 105, 

011302 (2010), 1003.2635. 
[38] C. Germani and A. Kehagias, JCAP 1005, 019 (2010), 

1003.4285. 

[39] J. Garcia-Bellido, D. G. Figueroa, and J. Rubio, Phys. 

Rev. D79, 063531 (2009), 0812.4624. 
[40] F. Bezrukov, D. Gorbunov, and M. Shaposhnikov, JCAP 

0906, 029 (2009), 0812.3622. 
[41] A. Zee, Phys. Rev. Lett. 42, 417 (1979). 
[42] S. Weinberg, Gravitation and Cosmology (Wiley, NY, 

1972). 

[43] N. D. Birrell and P. C. W. Davies, Quantum Fields in 
Curved Space (Cambridge University Press, Cambridge, 
UK, 1982). 

[44] K.-i. Maeda, Phys. Rev. D39, 3159 (1989). 

[45] M. E. Peskin and D. V. Schroeder, An Introduction to 

quantum field theory (Perseus Books Publishing, L.L.C., 

1995). 

[46] J. R. Ellis, M. K. Gaillard, and D. V. Nanopoulos, Nucl. 

Phys. B106, 292 (1976). 
[47] A. I. Vainshtein, M. B. Voloshin, V. I. Zakharov, and 

M. A. Shifman, Sov. J. Nucl. Phys. 30, 711 (1979). 
[48] J. Steinberger, Phys. Rev. 76, 1180 (1949). 
[49] L. Resnick, M. K. Sundaresan, and P. J. S. Watson, Phys. 

Rev. D8, 172 (1973). 
[50] T. G. Rizzo, Phys. Rev. D22, 178 (1980). 



